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Abstract: We study the motion of a quantum charged particle, constrained on the surface of a cyhnder, 
in the presence of a radial magnetic field. When the spin of the particle is neglected, the system essentially 
reduces to an infinite family of simple harmonic oscillators, equally spaced along the axis of the cylinder. 
Interestingly enough, it can be used as a quantum Fourier transformer, with convenient visual output. When 
the spin 1 /2 of the particle is taken into account, a non-conventional perturbative analysis results in a recursive 
closed form for the corrections to the energy and the wavefunction, for all eigenstates, to all orders in the 
magnetic moment of the particle. A simple two-state system is also presented, the time evolution of which 
involves an approximate precession of the spin perpendicularly to the magnetic field. A number of plots 
highlight the findings while several three-dimensional animations have been made available on the web. 
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1 Introduction 

The quantum mechanical description of the motion of charged particles in a magnetic field is a classic appli- 
cation of the theory, having given rise to some of its most striking results. Among them, the seminal analysis 
by Dirac 1^1, of the motion in the field of a magnetic monopole, continues to inspire decades after its incep- 
tion, and motivates the study of similar quantum systems that share the characteristic of providing insights 
into the fundamentals without too much distraction by analytical complexity. Such systems are invaluable 
pedagogically, as they furnish a manageable, yet captivating testing ground of the fundamentals of the theory. 

The problem of the motion of a non-relativistic quantum particle in a plane, in the presence of a perpen- 
dicular homogeneous magnetic field is presented in several textbooks (see, e.g., |13| ) — nevertheless, it seems 
to be the only standard example of this type available. The main purpose of this paper is to draw attention 
to the fact that the analogous problem for the cylinder is also manageable, even when augmented to include 
a spin 1/2. In this latter case, we also show how the use of the creation and anihilation operator machinery 
greatly simplifies the perturbative analysis of the problem, in comparison to the standard textbook procedure. 

Despite the simplicity of the problem and it being an obvious variation on the monopole theme, we have 
not been able to find a treatment in the literature. The motion of a spin-1/2 particle in the field of a magnetic 
monopole has been studied in detail, both in the non-relativistic [7| HI I15L and relativistic f2] cases. 
Symmetry aspects of the problem have also been considered extensively (see, e.g., ^3)) with the discovery 
of an underlying supersymmetry among the most notable results ^ El EI]- On the other hand, quantum 
spinless particles moving on curves or surfaces have been extensively studied (see, e.g., [31 I16L I14| and 
references therein) with a general discussion of the effects of a vector potential given in . It is our hope that 
the use of the above simple system will enhance the exposition of this fascinating part of the theory. It should 
also be of interest in practical applications, such as constrained quantum mechanics and carbon nanotube 
physics. 

Consider a classical charged particle, constrained to move on the surface of an infinite cylinder, in the 
presence of a radial magnetic field, 

B{r) = Bo-p, (1) 
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where a is the radius of the cyUnder and Bq is the field strength on its surface-*^. The equations of motion for 
the particle are 

mv^ = -qBo , mVif, = qB^Vz , (2) 

where to, q are the mass and charge of the particle respectively and {p, (f), z} is a right-handed basis. The 
solutions to Q are two simultaneous oscillations: the momentum pz of the particle oscillates like, say, cos{ujt) 
(with Lu = qBo/mc) while its angular momentum along the z-axis oscillates like sin(wi). Thus, the particle's 
kinetic energy oscillates between a linear and a rotational form, becoming, for example, purely rotational at 
the turning points of the oscillation along z. 

We study, in this paper, the quantum mechanical version of the above problem, adding, at a later stage, 
a spin-1/2 to the particle. The treatment of the spinless case, contained in Sect.|21 is exact — the problem 
separates and reduces to an infinite collection of harmonic oscillators along z. We find, nevertheless, the 
resulting quantum system particularly rich and with surprising properties — it functions, for example, as a 
quantum Fourier transformer with convenient visual output (see Sec. I2.2)l . The addition of spin is treated 
perturbatively in Sect. 13 with a non-conventional method that greatly simplifies the calculations. We are 
able to give recursion relations for the corrections to the wavefunctions and the energy to all orders, for 
all unperturbed eigenstates, and apply the results to compute second-order corrections to the ground state. 
Several plots highlight the findings. We also make available on the web several three-dimensional color 
animations of the time evolution of the wavefunction, with or without spin, and corresponding to various initial 
conditions. An appendix shows how the standard perturbation theory treatment of the problem reproduces, 
albeit laboriously, our first order results. 



2 The Spinless Case 
2.1 The spectrum 

The magnetic field of Eq. (QJ can be obtained, in the vicinity of the surface of the cylinder, from the vector 
potential 

A{r)=-Bo-z$. (3) 
P 

The Hamiltonian for a quantum spinless particle constrained to move on the surface of the cylinder is given 

by 

Zm c 



2m'a'^ ^' 2toc2 



i^&■i + ^t) + ^z^-^^z^,. (4) 



The wavefunction *(0,z) = -i=e^^'^Z(z) is an eigenfunction of , with eigenvalue E, provided Z(z) satisfies 
(primes denote differentiation w.r.t. z) 

-^Z"{z) + \mu?{z + ihfZ{z) = EZ{z) , (5) 
2to 2 

where we have set 

, b= , (6) 

mc qBoa 

and, in what follows, we take h = m — lu ~ 1. This is the equation for a simple harmonic oscillator (SHO), 
centered at z = —£b. We conclude that, for each integer value of £, one obtains a copy of the usual SHO 



'^Such a radial field can be thought to be produced by a homogeneous linear magnetic charge density, or, more realistically, 
in the exterior of a solenoid of radius R, R < a, placed along the axis of the cylinder and carrying surface current density 
J = -^^z4>. 
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spectrum, centered at z = —ib, i.e., the eigenfunctions and eigenvalues of H are given by 

(0, z\n,e)^ NMz + ^6)e-(-+^'')'/2 1 ^ e.^^ = E„ = n + ^ , iV„ = (2«n!V^)-i , 

(7) 

where \n, £) denote the corresponding eigenkets (rt = 0, 1, 2, . . . ; £ S Z) and Hn{z) are the Hermitc polyno- 
mials. 

2.2 A quantum Fourier transformer 

Suppose that the wavefunction of the particle, at i = 0, is given by^ 

^e{4>, z, i = 0) = (0, z I 0, 0, e) = Noe-'"/^^e'"^ . (8) 

If ^ = 0, we have one of the infinitely many ground states of the system and the time evolution is by a phase 
factor. Consider now the case i ^ 0. Then the z-part "sees" a quadratic potential centered at z = —ib but 
the initial wavefunction is a gaussian centered at the origin. This is a coherent state and its time evolution is 
an oscillation around z = —ib, with the frequency w = 1 of the oscillator, 

Z, t) = iVoe-^*/' ^~^(z+a)lbsint^^{z+eb{l-cost)f/2 1 _ (9) 

V 27r 



Notice that a physically irrelevant global phase factor e'^^^^/^^'"^* j-^g^g i-,ggu omitted from the above expression. 
We may now exploit linearity to write down the time evolution of a gaussian (in z), centered at the origin, 
with arbitrary 0-dependence, 

^{<l>, z, t = 0) = iVoe-^'/2/(0) • (10) 



Writing 



we obtain, 

oo 

*(</), z,i)= J2 fi'^ii<l>, (12) 

l=-oo 

i.e., each Fourier mode of f{(f)) gives rise to a gaussian in z, oscillating like a coherent state around z = —£b 
with frequency lo = 1. Taking 6 3> 1, so that the various gaussians separate after half a period, converts 
the system to a quantum Fourier transformer with convenient visual output: looking at the wavefunction at 
time t — TT {Si half-period), one sees the above gaussians at the (second) turning point of their oscillation, at 
z = —2£b, with their amplitudes proportional to the Fourier amplitudes fg. In Fig.^ we plot several frames 
of the time evolution of '5, when /((/>) = cos0 — the last frame, at t = tt, clearly displays the Fourier content 
of /. Fig. 121 corresponds to the initial wavefunction *((/), z,t^O) ^ e~^'/2(l + e"'"^ -f- l.Se^*^'^ -t- £"'30), 



3 The Spin-1/2 Case 
3.1 Separation of variables 

For a spin-1/2 particle, the wavefunction has two components, ^'+(0, z), z), which we arrange in a 

column vector. The spin interacts with the magnetic field via Hint ~ —AS* • B, which, for the field given in ^ 

^We use the notation (0, z\n, £, m) = N„Hn{z + £fe)e~(^+^'') /2_^_g»"n^> (jio\^ to be confused, hopefully, with the standard 

spherical symmetry notation) — these wavefunctions are eigenfunctions of H only when m = l,\a which case they will be denoted 
by |n, £), as above. 
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Figure 1: Time evolution of the initial wavefunction '3/(0, z, t — 0) — 7Voe~^^/^ cos 0. The modulus 
of 'J is indicated by the radial distance of the surface from that of the cylinder while its phase is color-coded, with 
1, i, —1, — i corresponding to red, green, blue, purple (several animations in color, including the above, can be 
seen at http: //www.nuclecu.unam.mx/^chryss). The time t is equal to zero at the top left and increases to 
the right and downwards, reaching t = tt (half a period) at the bottom right. 



becomes 



1 







,-i<t> 





In order to achieve separation of variables now, we need to take 

1 



b, z) 



12-K 



The resuhing equation for the Z's is 



27r 



(13) 



(14) 



HiZ+{z) + €Z-{z) = EZ+{z) 
He+iZ^{z) + eZ+{z) = EZ_{z) , 



where Hi is a SHO hamiltonian centered at z = —ib 

1 



(15) 



(16) 
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Figure 2: Time evolution of the initial wavefunction *(</>, z, t = 0) ~ e~^^/^(l + e~*'^ + 1.5e~^^'^ + 
g— 13</)^_ Conventions are as in Fig.U Notice how the amplitudes of the gaussians in the last frame correspond 
to the Fourier components of the initial wavefunction. 



and e = —hXBQ/2. We see that the problem reduces to that of two SHO's, a distance b apart, coupled by 
the e terms in (|15ll . Our task is to solve p5(l perturbatively in e. Once the solutions are known, to a certain 
order in e, we can form the spinor 



^) = ,r !' - , (17) 



from which the probability density and spin direction can be extracted as 



\Z-\ „ . „ 



p= Vl^+P + I^-P, a = 2arctan^, /3 ^ im (log — ) , (18) 

|Z+| v+ 

where the spin direction n is given by n = (sin a cos /3, sin a sin /3, cos a), in Cartesian coordinates. The 
solutions of (|15|l have no relative (complex) phase and can be taken real. Then, on the (p — plane, the 
spin lies in the x — z plane. The extra e"^ factor in guarantees that when we change our position on the 
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cylinder by 0, the spin also rotates by the same angle and its direction is therefore obtained by revolution of 
the (f) = configuration, in other words, P = (f) in (|18|l . These remarks are of course valid only for the energy 
eigenstates — the time evolution of general states results in the spin pointing outside of the radial plane as 
well, even if they start within it. 



3.2 Perturbative solution 

One may treat the system of the two coupled differential equations in (|15|) by standard perturbation theory 
methods — we give the first-order analysis along these lines in the appendix. It is instructive though, as 
well as much more efficient, to exploit the SHO machinery of raising and lowering operators. We begin 
by transforming H15|) into a single differential-difference equation. Indeed, it is clear from the symmetry of 
these equations that the solutions can be taken to satisfy Z-{z) = ±Z^(—z — b) — we will refer to the two 
possibilities as symmetric and antisymmetric respectively. Taking t — Q and restricting to the symmetric case, 
the first of p5|l becomes 

HoZ,+ {z) + eZ,+ {^z - 6) = EZ,+ {z) . (19) 
We now write H{) = a^a + ^ and introduce the ket whose wavefunction is Zs+(2:), 

oo 

Z,+ {z)^{z\Z,+) , |Z,+) = ^c„(at)"|0)o = /(at)|0)o, (20) 

where |0)o is the ground state for t equal to zero {i.e., centered at the origin) and /(a^) is defined by the last 
equation — the idea is that any ket can be obtained by some function of applied to the ground state. The 
wavefunction Zs+{—z — b) is obtained from Zs+{z) by first reflecting around the origin and then effecting the 
translation z i-^ z + b. Reflecting around the origin an eigenfunction of Hq introduces a sign given by the 
parity of the state, which shows that the reflected state is produced by f{—a'^) applied to the ground state. 
The reflected and translated state then is given by 

T_b/(-at)|0)o = e-^'/4e-75°'e7fV(-a^)|0)o 

= e-^y^e-^'^'fi-a^-^)\0)o, (21) 

which brings H19|l into the form 

(a+a + i)/(at)|0)o + ee^"' /'e~^''' f{-a^ ~ -|)|0)o = Ef{a^)\0)o . (22) 

Noting that the a — commutation relations are identical to the dx — x ones, we may infer a differential - 
difference equation for the function /(cc), 

^n^) + \f{^) + ee-^'/^e-^yi-x " ^) = Ef{x) . (23) 

We now specify to the case where the unperturbed state is the n-th excited state of the SHO. The perturbed 
state will be denoted by /n(a^)|0)o, with energy En, where 



oo oo 



fc=0 fc=0 

Notice that 



f(o)(:c)^^x-, E(^)=n+l. (25) 
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Substituting these expansions in H23|l we obtain 



+ (26) 

m—l 

where we separated the m = 0, k terms in the sum on the r.h.s. and used H25|l . Notice that the r.h.s. above 
only contains with m < k, so H2fi|l can be used recursively to determine any fn^\x). 

The requirement that the perturbed eigenket be normalized implies that the corrections, order by order in 
e, have to be orthogonal to the unperturbed eigenket |n)o. This in turn implies that, for k > 0, the coefficient 
of a;" in fii'\x) must vanish. Then so does the coefficient of x" in xdxf^\x). Using this information, we 
can extract the coefficient of on both sides of — the resulting equation fixes recursively the energy 
corrections eI^\ 

= (") ( " ^r'mit~n-^) . (27) 

This is an appropriate point to comment on the antisymmetric solutions. The difference in this case is that 
the e term in (|23|l appears with a minus sign. Since this is the only place where e appears explicitly, we 
conclude that one gets the antisymmetric solutions from the symmetric ones by the substitution e — e. As 
we will see later on, symmetric solutions have their spin parallel, more or less, with the magnetic field while 
antisymmetric ones have it antiparallel. Since e is proportional to the magnetic moment of the particle, the 
above statement about the relation between the two kinds of solutions essentially says that the symmetric 
solution for a particle coincides with the antisymmetric solution for the same particle but with the opposite 
magnetic moment. 

3.3 Corrections to the ground state 

For n = 0, Eqs. (|?7ll simplify considerably, 

^ m—l 

^(.) ^ ^-.V4;('=-i)(__L). (29) 

The unitarity argument given above implies in this case that /o'^^(O) = 0, for all k greater than zero — this 
fixes the lower integration limit in the solution of (|28|) equal to zero. The change of variable xj, xsk and 
the substitution of (|29|l finally give 

(30) 

Notice that the (apparent) pole of the integrand at Sfc = cancels out. 

We look in some detail now at the wavefunctions and resulting spin configurations, including up to 
quadratic corrections. For the first three f'^^\ Eqs. H25|) . H3U|I give 

/r(^) = 1 (31) 

i'\x) - e-^/* r^(l-e-^-) (32) 

Jo Jo ■S2S1 L J 
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while for the corresponding energy corrections we get 



E^^^ = I , eI,'^ = e-^'l' , Ef = 6-^^/2 ^ (l - e^-) . (34) 

JO ^1 ^ ' 



Applying the above /'s to the ground state and projecting onto the position eigenket \z) we find the wave- 
functions 

(35) 

^ dsi 



!i (e-.V2 „ gfc=.?/4g-(.+bs,)V2) (36) 
ds2dsi 



s 



■S2S 



r{0 



^ gf,2si/2+62s^(l-si)V4g-(z+6s2(l-si))V2 _ ^b^s^s? /4g- (z+bs2Si ) ^2 | ^ (37) 

In terms of these, the symmetric solution for the spinor has components (we ommit an overall normalization 
factor) 

Z,+ (z) = zf (z) + eZ«(z) + e2z|2)(z), Z^^) ^ Z,^{-z - h) (38) 

and energy 

E,, = Ef ^eE^^'> ^e^E^ , (39) 
while the antisymmetric solution is given by 

Z^^{z) = -Zf{z) + eZ%\z)-^Z%\z), Z,^{z) = -Z,+ {-z-b) (40) 

with energy 

E^„ = - .41) + .^4^) . (41) 

Plots of Z±, for both cases, as well as the corresponding spin configurations, are given in Figs. 01 E] Notice 
that the effect of the perturbation is small, despite a rather large value of e. This can be traced to the fact 
that the parameter values used give rise to a small overlap of two neighboring gaussians. As a result, the 
range of validity of our perturbative results is considerable larger that the standard e ^ 1. 



3.4 A two-state system 

Consider a state that, at t = 0, is the sum of the symmetric and antisymmetric n — states found above, 
for, say, £ — 0. The spin-down component of this state is of order 1 while the spin-up component is of order 
e. Roughly speaking, the particle is localized at z = —b and its spin points along the negative z-axis. Then 
the standard two-state system analysis shows that the amplitudes to be in the spin-up and spin-down states 
at later times behave like —i sin(rii) and cos{nt) respectively, where 's expansion in powers of e starts with 
eE^^^/2. The spin precesses in the tangent plane to the cylinder {i.e., perpendicularly to the magnetic field) 
while the particle oscillates from — & to zero and back. 



4 Concluding Remarks 

We have studied the problem of the motion of a spin-1/2 particle on a cylinder, in the presence of a radial 
magnetic field. A non-standard perturbative analysis, applicable to any perturbation of the harmonic ascil- 
lator, led to a recursion relation for the wavefunction and energy corrections, Eqs. (|26|l and (|27|l respectively, 
with explicit results for the ground state in Eqs. (|34|l - H37|l and Figures Ol El It is worth emphasizing that 
the radial magnetic field of the problem can be approximated in the laboratory, as pointed out in the first 
footnote. 
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z = -b z = 



Figure 3: The spin components Zs-\-, Zs— and the corresponding spin configuration. The dotted curves 
give the zeroth-order result, i.e., two gaussians, centered at z = and z = —b = —2. The dashed curves include 
corrections up to first order while the solid ones up to second (e ~ 0.5 — a rather large value was used to make 
the effect visible). The integrals in Eqs. p6(l . p7(l have been evaluated numerically. The spin configuration shown 
includes quadratic corrections. Due to the symmetry, the spin points always along p (upwards in the figure) at 
z — —6/2 = —1, while it tends to ±z as z tends to ±oo. The perturbation causes a zero in each component (for 
finite z) — the spin crosses the z axis at those points. As a result, the two gaussians are pushed apart while the 
width of the central region, where the spin points up, is reduced. 



We end with a comment on the form of the unperturbed hamihonian used, Eq. When deahng with 
the motion of a quantum particle on a surface, one can use a 3-D Laplacian in the hamiltonian and constrain 
the motion of the particle on the surface using a steep confining potential in the radial direction. It is well 
known that, in this approach, which seems to be the one appropriate for practical applications, there is an 
induced potential for the motion along the surface, proportional to the square of the difference between the 
two principal curvatures of the surface (see, e.g.., p^ ll4j and references therein). In our case, this is a constant 
which only shifts the energy eigenvalues. Nevertheless, an obvious extension of our problem here would be the 
study of the motion on the surface of a slightly curved cylinder, in which case the above mentioned induced 
potential would have to be taken into account. 



A First Order Corrections to the Ground State: the Standard 
Treatment 

We outline here the standard first order perturbative analysis of the problem, deriving the corrections to 
the ground state. Given that the zeroth order spectrum is degenerate, we need to first diagonalize the 
interaction hamiltonian in each degenerate subspace. One easily sees that the interaction only connects the 
pairs of eigenstates |n, £, +), |n, £ +1, — ). The appropriate zeroth order basis is given by the symmetric and 
antisymmetric linear combinations 

In?) ^ K) = ^{\n,£,+) + \n,e+l,-)) 

K> ^ K> - ^{-\n,£,+) + \n,£+l,-)). (42) 

Notice that we use the label £ for states that are equally localized at —£b and ~{£+l)b. The reason for renaming 
the states with numerical superscripts, instead of letters, will become apparent below. The expectation value 
of iJint in these states reproduces our result (|34|l for the first order correction to the energy. It is interesting 
to see how the first order correction to the wavefunction, conventionally given by an infinite sum, is brought 
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z = -h z = 



Figure 4: The spin components Za+, Za_ and the corresponding spin configuration. Notation and 
parameter values are as in Fig. |31 Due to the antisymmetry, the spin points always along — p (downwards in the 
figure) at z = —5/2 — —1, while it tends to ±z as z tends to ±oo. The perturbation pushes the two gaussians 
together while the width of the region where the spin points downwards is increased. 



into the closed form . The matrix elements of Hint in the above basis are 

(n°|iJint|m° ) = e {ni\me+i) Su'Shrh 
{n1\Hint\m},) = e {ne\me+i) 6u'6ii,rh+i 

{nl\Hint\ml,) = -e{ni\mi+i)Sii'Snm, (43) 

where {n{ \ m^+i) is the overlap between SHO eigenstates |m), at a distance b apart and n is the parity 
of n. Specifying to the symmetric ground state and taking £ = 0, we find the first order correction 

100^(1) ^gMMlfc^^). (44) 
fe=i 

We see that only states with £ = contribute. Furthermore, when k is even, only the symmetric state 
contributes while for k odd, only the antisymmetric one does. Using the fact that 

(Oo|fci) = -^(--^)'^e-^, (45) 



we find for the spin-up component of the correction 



which implies 



= (|0i) - |0o)) . (47) 

Integrating back w.r.t. b and changing integration variable we recover our earlier result (|32|l (similarly for the 
spin-down component). 
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